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Analysis of parametric oscillatory instability in Fabry-Perot cavity with Gauss and 

Laguerre- Gauss main mode profile 
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We calculate the parametric instabilities in Fabry-Perot cavities of Advanced VIRGO and LIGO 
interferometers with different main mode profiles. All unstable combinations of elastic and Stokes 
modes both for the case with TEMoo and LG33 as a carriers are deduced. 

PACS codes: 42.50.Wk, 42.65.Es, 43.40.+S, 46.40.-f 

Keywords: parametric oscillatory instability, LIGO and VIRGO interferometers 

I. INTRODUCTION 

Sensitivity increase of gravitational wave detectors like Advanced VIRGO and LIGO is expected to be obtained by 
rising the value of power stored in the Fabry-Perot(FP) resonator optical mode. However, high values of circulating power 
may be a source of the nonlinear effects which will prevent from reaching the projected sensitivity. It is appropriate to 
remind that nonlinear coupling of elastic and light waves in continuous media produces Mandelstam-Brillouin scattering. 
It is a classical parametric effect; it is often explained in terms of quantum physics: one quantum Tiwo of main optical 
wave transforms into two, i.e., Tuui in the additional optical wave (Stokes wave: cui < cuq) and ?uu m in the elastic wave 
so that Wo — <JJ] + cu m . The irradiation into the anti-Stokes wave is also possible (u)i a = cuq + cu m ), however, in this 
case the part of energy is taken from the elastic wave. The physical mechanism of this coupling is the dependence of 
refractive index on density which is modulated by elastic waves. If the main wave power is large enough the stimulated 
scattering will take place, the amplitudes of elastic and Stokes waves will increase substantially. 

In gravitational wave detectors elastic oscillations in FP resonator mirrors will interact with optical ones being coupled 
parametrically due to the boundary conditions on the one hand, and due to the ponderomotive force on the other hand. 
Two optical modes may play roles of the main and Stokes waves. High quality factors of these modes and of the elastic one 
will increase the effectiveness of the interaction between them and may give birth to the parametric oscillatory instability 
which is similar to stimulated Mandelstam-Brillouin effect. This undesirable effect of parametric oscillatory instability 
may create a specific upper limit for the value of power W c circulating inside the cavity Q. 

It is interesting that the effect of parametric instability is important not only for large scale gravitational-wave in- 



terferometers. Recently, the instability produced by optical rigidity was observed in experiment 



El- K. 



Vahala with 



collaborators has also observed it in micro scale whispering gallery optical resonators [3|, [4|]. Zhao et al.[5j have shown 
in experiment, using an 80 m Fabry-Perot cavity, that parametric instability effect can indeed occur. The experimental 
results are in good agreement with theoretical predictions 1 1 . 



The condition of parametric instability for Fabry-Perot cavity may be written in simple form[l(: 
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where c is a speed of light. L is a distance between the FP mirrors, m is a mirror's mass, y m , yi are the relaxation rates of 
clastic and Stokes modes correspondingly, W is a power circulating inside the cavity, Ai is an overlapping factor of clastic 
and optical modes and A = cuo — u^i — w m is a detuning value. Here Ao and Ai are the functions of the distributions 
over the mirror surface of the optical fields in the main and Stokes optical modes correspondingly, vector u is the spatial 
vector of displacements in elastic mode, u z is the component of u, normal to the mirror's surface, J dr± corresponds to 
the integration over the mirror surface and J" dV - over the mirror volume V. 

As we can see from formula ([1]) the parametric instability is a threshold effect and it takes place if optical power W 
in main mode of Fabry-Perot cavity is bigger than the threshold power W c . The value W c depends on the detuning A 
between optical and elastic modes and overlapping factor Ai , in particular, the W c has minimum if detuning is less than 
relaxation rate of Stokes mode A yi and overlapping factor is large enough. 




FIG. 1: Schematic structure of optical Laguerre-Gauss modes in the FP cavity. The modes of the main frequencies sequence are 
shown by higher peaks. It is shown that Stokes mode with frequency u)i may not have suitable anti-Stokes mode (it is denoted by 
question-mark) . 

Analyzing this equation D'Ambrosio and Kclls 6] pointed out that the presence of the anti-Stokes mode can considerably 
depress or even exclude parametric instability. Indeed, the condition of parametric oscillatory instability in the presence 
of anti-Stokes mode with frequency cui a = cuq + cu m has the following form[7[: 



where y-\ a is a relaxation rate of anti-Stokes mode and Ai a = Wi a — cuo — cu m is a detuning value. 

For example, let the main, Stokes and anti-Stokes modes be equidistant and belong to the main frequency sequence 
cui = tc(K — 1)c/L, o>o = 7tKc/L, (Via — 7t(K + 1 )c/L(K is an integer). In this case A = Ai a , the main, Stokes and 
anti-Stokes modes have the same Gaussian distribution over the cross section and hence the same overlapping factors: 
Ai = Ai Q . It means that the second term in the right part of equation is larger than first term, the positive damping 
introduced into elastic mode by the anti-Stokes mode is greater than negative damping due to the Stokes mode, hence 
the parametric instability is impossible. This case has been analyzed in details in 6|. 

For the case when the Stokes and anti-Stokes modes do not belong to the main sequence (non-zero radial N and 
azimuth M numbers) the frequencies of the suitable Stokes and anti-Stokes modes are not equidistant from the main 
mode (A ^ Ai Q ) and have different spatial distributions (A] ^ Aj Q ). It takes place when main mode is TEMqo and Fig. 
Q] illustrates it. For the shown Stokes mode (left to the main mode) there is no suitable anti-Stokes mode (it should be 
located right to the main one). The probability that appropriate anti-Stokes mode exists is extremely small and we see 
that in this case full depression of parametric instability docs not take place. However, when main mode has non-zero 
radial and azimuth numbers(high order Laguerre-Gauss modes LG pn ) the possibility of parametric instability depression 
by existence of suitable anti-Stokes mode will increase. 

Parametric instability is a serious problem for advanced gravitational wave detectors and we need to know this "enemy" 
in detail to avoid it. In particular, we need detailed information about pairs of Stokes modes and elastic modes which 
may be possible candidates for parametric instability. It is known that Stokes and anti-Stokes modes may be analytically 
calculated for Gaussian beams. In contrast, elastic modes may be calculated numerically(some method of analytical 
clastic modes calculation has been proposed in [§|). Hence, the accuracy of parametric instability forecast (and success 
of methods to prevent it) directly depends on how accurately we can calculate normal (eigen) frequencies and spatial 
distributions of elastic modes. Using ANSYS® code, for example, we can obtain the accuracy of elastic modes calculations 



of about 0.5% or poorer [9 
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In [131 141 the importance of accurate numerical calculation of clastic modes in the mirrors of Advanced LIGO was 

n 

discussed to enable precise predictions of the problem of parametric oscillatory instability. In [13[ it has been proposed 
accuracy estimations through use of analytical solutions based on Chree-Lamb modes. Small deviations from cylindrical 
test mass shape may produce splitting of non-axial symmetric elastic modes into doublets. This splitting may increase 
the possibility of parametric oscillatory instability too 14 1. 



It is known the important goal in gravitational wave detectors is to minimize fundamental and technical noise contri- 
bution. For example, progressive way to lower the thermal noise is to change the mode shape of the laser beam inside 
the interferometer. It is worth noting that different shapes of beam have been proposed for reducing thermal noise such 



as mesa beams 
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17| . conical modes 



18[ and high order Laguerre-Gauss modes 



techniques for the generation of hi gh o rder LG modes using holograms 



with high conversion efficiency 
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19(. At present, there are some 



21| , gratings 22| and mode transformers 
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26j using numerical interferometer simulations the comparison of behaviour of 
the LG33 mode with the fundamental mode TEMoo are discussed. The LG33 mode performs similar if not even better 
than commonly used TEM.00 for all considered aspects of intcrfcromctric sensing. 

In this Letter we estimate the possibility of parametric instabilities in Fabry-Pcrot cavity of Advanced VIRGO and 



26| . In section [TT| we provide 



LIGO interferometers both for Gauss TEMoo and Laguerre-Gauss LG33 modes as a carriers 
estimations of PI in Fabry-Perot cavity of Advanced VIRGO and LIGO interferometers and discuss our results, in section 
IIIII we make some useful conclusions for future research. 



II. RESULTS 

In order to predict the unstable combinations of Stokes and elastic modes we have to take into account the additional 
azimuth numbers condition of parametric instability. For example, for TEMoo as a carrier (when the elastic modes have 
azimuth dependence - e 1 ™* and Stokes modes have - e ln *) the non-zero overlap factor will be if m = n(see the formula 
([2]) for Ai ). In turn, if we use LG33 as a carrier the non-zero overlap factor will be if |m ± n| = 3. It is worth noting 
that it is correct only if cylinder center coincides with laser spot center and below we consider this particular case (in 
opposite case the overlap factor has to depend on distance between the center of mirror and center of main optical mode 
distribution over mirror surface). 

For Fabry-Perot cavity we use the parameters of Advanced VIRGO and LIGO interferometers shown in Table|T]where T 

TABLE I: Parameters of Advanced LIGO and VIRGO interferometers 



VIRGO 


LIGO 


W = 0.76 x 10 6 W 


W = 0.83 x 10 6 W 


L = 3000m 


L = 4000m 


T = 7 x 10~ 3 


T = 14 x 10~ 3 


A = 1064nm 


A = 1064nm 


m = 40kg 


m = 40kg 


r = 0.17m 


r = 0.17m 


H = 0.2m 


H = 0.2m 


y m = 6 x 10 _3 s _1 


y m =6 x 10~ 3 s~' 


y = cT/4L = 175s" 1 


y = cT/4L = 262.5s" 1 



are the transmittances of input FP mirrors, A is a wave length of input laser and r, H are radius and height of cylindrical 
fused silica mirror. Note that the condition of parametric instability in Michelson interferometer with Fabry-Pcrot 
cavities in the arms with additional power and signal recycling mirrors (configurations of Advanced LIGO and VIRGO 
interferometers) and in alone Fabry-Perot cavity can substantially differ - see discussion in Conclusions. 

The calculation of parametric instabilities is very sensitive to small changes in the model parameters and gives only 
statistical result for number of unstable modes. In particular, uncertainties in the mirror radius of curvature and small 
variations in materials and temperature are sufficient to move the optical and mechanical resonances of the cavities. We 
estimate the number of unstable combinations of elastic and Stokes modes both for the case with TEMoo as a carrier(TEMoo 
case) and - LG33 as a carrier(LG33 case). The elastic modes were calculated numerically using COMSOL® code on triangle 
mesh with about 40000 meshing elements. We restrict our estimates by elastic frequency range up to 2n x 40000s 
because, on the one hand, the more the elastic frequency the less the parametric gain. On the other hand, high elastic 
mode distributions arc difficult to determine exactly among other elastic modes due to degradation of numerical accuracy 
of COMSOL® code. 

For VIRGO configuration we use Fabry-Perot cavities with laser spot radii on mirrors surfaces w = 6.47 cm and 



w = 3.94 cm in TEMoo and LG33 cases correspondingly [26|. The waist radii are situated in the middle of the cavities 
and are equal to wq ~ 0.79cm and Wq ~ 1 .38cm for these cases. Level of diffractional losses in clipping approximation is 
hold to be l c up — 1 ppm. For LIGO configuration we use Fabry-Perot cavities with laser spot radii on mirrors surfaces 



27j. The waist radii are situated in the middle of 



w — 6 cm and w = 3.94 cm in TEMoo and LG33 cases correspondingly : 
the cavities and are equal to wo ~ 1.15cm and wo ~ 1.99cm for these cases, and level of diffractional losses in clipping 
approximation is hold to be l c u p = 1 ppm or less. 

We calculate the values of overlap factors for all combinations of elastic and optical modes(up to 9 th order) but, on 
the other hand, we use azimuth numbers condition in each case. It is also worth noting that the analysis of parametric 
instability in each case has been realized with taking into account the appropriate anti-Stokes modes which may suppress 
the parametric instabilities. We also took into account the diffractional losses of Stokes and anti-Stokes modes in clipping 
approximation for estimations of optical modes relaxation rates as it has been done in [7|. 

In Table HT1 all unstable combinations of elastic and Stokes modes and parametric gains R for LG33 case in VIRGO 
interferometer are shown. There is only one unstable combination and there is not any anti-Stokes mode that can suppress 
the parametric instability. In figurc[5]we can see the z-component displacement vector distribution for this unstable clastic 
mode. In the case of TEMoo carrier we did not find any unstable modes in given frequency range due to small overlap 
factors and large detuning values A. 

In Table IIIII all unstable combinations of elastic and Stokes modes and parametric gains R for LG33 case in LIGO 



FIG. 2: Displacement vector component distribution u z of elastic mode with frequency 38577Hz and azimuth index m = 1 . The 
numerical calculations of cylinder have been made on triangle mesh with about 40000 meshing elements. 

interferometer are shown. There are three unstable combinations. In figure [3] we present the z-component displacement 
vector distributions for these elastic modes. In the case of TEMoo carrier we found only one unstable mode in given 
frequency range which is shown in Tabic HVl 

We can see that possibility of parametric instability in given elastic frequency range for LG33 case will be slightly larger 
than for the TEMoo case both in Fabry-Perot cavities of VIRGO and LIGO interferometers. 

TABLE II: Unstable combination of elastic and Stokes optical modes in FP cavity of VIRGO configuration with LG33 as a carrier. 
In the case of R1M.00 carrier we did not find any unstable modes in given frequency range due to small overlap factors and large 
detuning values A. 



w m /2n, Hz 


m 


Stokes optical mode 


R 


38577 


1 


LG32 


2.5 



TABLE III: Unstable combinations of elastic and Stokes optical modes in FP cavity of LIGO configuration with LG33 as a carrier. 



cu m /27t, Hz 


m 


Stokes optical mode 


R 


36524 


6 


LG 03 


1.5 


37566 


4 


LG 17 


9.3 


37566 


4 


LG 41 


10.8 
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FIG. 3: Displacement vector component distributions u z of elastic modes with frequencies 23048Hz (azimuth index m = 1), 
36524Hz (azimuth index m = 6) and 37566Hz (azimuth index m = 4) correspondingly. The numerical calculations of cylinder have 
been made on triangle mesh with about 40000 meshing elements. 

TABLE IV: Unstable combination of elastic and Stokes optical modes in FP cavity of LIGO configuration with TEMoo as a carrier. 



w m /2n, Hz 


m 


Stokes optical mode 


R 


23048 


1 


LG21 


1.7 



III. CONCLUSIONS 

In this Letter we deduced the analysis of parametric oscillatory instability in FP cavities of gravitational wave detectors 
like Advanced LIGO and VIRGO in the model with cylindrical mirrors without flats and suspension ears. The main goal 
of our estimations was to compare a difference in TEM00/LG33 behavior between Advanced LIGO and VIRGO Fabry- 
Pcrot cavities for parametric instabilities in such simple model. The number of obtained unstable modes for LG33 case 
is slightly larger than for TEMoo case in elastic modes frequency range up to 40kHz, however, this difference is not large 
enough to conclude that LG33 case is more dangerous than TEMoo case. No doubts, there is necessity to perform detailed 
analysis for full scale schemes of Advanced LIGO and VIRGO in LG33 ca se with wide elastic modes frequency range, 



very high accuracy of elastic modes calculations 



and mirror flats and ears. 
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At the end of the Letter we would like to summarize some dangerous points that must be taken into account at present 
and in future parametric instability research. 

It is worth noting that for combinations of modes suitable for parametric instability the overlapping factor Ai may 
be zero (for example, elastic mode and the Stokes mode can have different dependence on azimuth angle) — we did not 
mention such combinations in our estimations. However, it is important to take into account that only the elastic mode is 
attached to the mirror axis in contrast to the optical mode which can be shifted from the mirror axis due to non-perfect 
optical alignment. Hence, the overlapping factor has to depend on distance Z between the center of mirror and the center 
of the main optical mode distribution over the mirror surface. It means that Ai may be zero for Z = but non-zero for 
Z ^ 0. Therefore, the numerical analysis of the mode structure should evidently include the case when Z ^ 0. Note that 
there is a proposal to use special shift Z of the laser beam of about several centimeters from the mirror axis in order to 
decrease thermal suspension noise [28 1. 

Imperfections from the cylinder shape (such as flats and suspension ears), should cause splitting of elastic modes with 
azimuth index m > 1 into doublets and the difference between doublet frequencies may be large enough or greater as 
compared with relaxation rate of Stokes mode. Both the appearance of doublets and high density of elastic modes may 
increase the possibility of parametric instability in Advanced VIRGO and LIGO 



13 



14] 



The effect of parametric instability for power recycled interferometer may be larger than for the separate Fabry-Pcrot 
cavity because the Stokes mode emitted from the Fabry-Perot cavity throughout its input mirror is not lost irreversible 
but returns back due to power recycling mirror, therefore, its interaction is prolonged 7[. As it was demonstrated in 
29j the presence of power recycling cavity causes parametric gains R to be hugely amplified at A = and reduced by a 
factor of 2 in off-resonance intervals as compared to an interferometer without power recycling. The detailed analysis of 
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31 



On the other hand, the parametric 



parametric instability in signal recycled interferometer has been realized in 
instability in the interferometer with pow er recycling mirror may be depressed by variation of distance between power 
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33 
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Detail analysis of parametric instability in full scale Advanced LIGO 



34| and it will be useful to present the same analysis for LG33 case in future 



recycling mirror and beam splitter 
interferometer has been performed in 
researches. 

We think that parametric oscillatory instability effect can be overcome in laser gravitational detectors after developing 
strategies for its suppression using these investigations. 
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